Astronomy &: Astrophysics manuscript no. 19253 


© ESO 2013 


March 4, 2013 





Amplitudes of solar-like oscillations in red giant stars 
Evidence for non-adiabatic effects using CoRoT observations 

R. Samadi 1 , K. Belkacem 1 , M.-A. Dupret 2 , H.-G. Ludwig 3 * 4 , F. Baudin 5 , E. Caffau 3 ' 4 , M.-J. Goupil 1 , and C. 

Barb an 1 



(N 

o 

(N 



6 

> 

(N 
> 

00 

in 
o 

> 

X 

C3 



LESIA, CNRS UMR8109, Observatoire de Paris, Universite Pierre et Marie Curie, Universite Denis Diderot, Place 
Jules Janssen, 92195 Meudon Cedex, France 

Institut d'Astrophysique et de Geophysique de l'Universite de Liege, Allee du 6 Aout 17 - B 4000 Liege, Belgium 
Zentrum fur Astronomie der Universitat Heidelberg, Landessternwarte, Konigstuhl 12, D-69117 Heidelberg, Germany 
GEPI, CNRS, Observatoire de Paris, Universite Denis Diderot, Place Jules Janssen, 92195 Meudon Cedex, France 
Institut d'Astrophysique Spatiale, CNRS, Universite Paris XI, 91405 Orsay Cedex, France 



March 4, 2013 



ABSTRACT 



Context. A growing number of solar-like oscillations has been detected in red giant stars thanks to CoRoT and Kepler 
space-crafts. In the same way as for main-sequence stars, mode driving is attributed to turbulent convection in the 
uppermost convective layers of those stars. 

Aims. The seismic data gathered by CoRoT on red giant stars allow us to test mode driving theory in physical conditions 
different from main-sequence stars. 

Methods. Using a set of 3D hydrodynamical models representative of the upper layers of sub- and red giant stars, 
we computed the acoustic mode energy supply rate (CP max ). Assuming adiabatic pulsations and using global stellar 
models that assume that the surface stratification comes from the 3D hydrodynamical models, we computed the mode 
amplitude in terms of surface velocity. This was converted into intensity fluctuations using either a simplified adiabatic 
scaling relation or a non-adiabatic one. 

Results. From L and M (the luminosity and mass), the energy supply rate ^Pmax is found to scale as (L/M) 2 ' 6 for both 
main-sequence and red giant stars, extending previous results. The theoretical amplitudes in velocity under-estimate the 
Doppler velocity measurements obtained so far from the ground for red giant stars by about 30 %. In terms of intensity, 
the theoretical scaling law based on the adiabatic intensity-velocity scaling relation results in an under-estimation by 
a factor of about 2.5 with respect to the CoRoT seismic measurements. On the other hand, using the non-adiabatic 
intensity-velocity relation significantly reduces the discrepancy with the CoRoT data. The theoretical amplitudes remain 
40 % below, however, the CoRoT measurements. 

Conclusions. Our results show that scaling relations of mode amplitudes cannot be simply extended from main-sequence 
to red giant stars in terms of intensity on the basis of adiabatic relations because non-adiabatic effects for red giant 
stars are important and cannot be neglected. We discuss possible reasons for the remaining differences. 

Key words, convection - turbulence - atmosphere - Stars: oscillations - Stars: red giants 



1. Introduction 

Before CoRoT (launched in December 2006), solar- like os- 
cillations had been detected for a dozen of bright red gi- 
ant st ars either from the ground or from space with MOST 
(e.g., |Barban et al.|2007[ ). Thanks to CoRoT and Kepler, it 
is now possible to detect and measure solar-like oscillations 



in many more (several thousands) red 


giant stars (e.g., 


de 


Ridder et al.||2009| Huber et al.||2010 


Bedding et al.||20 


L0 


Kallinger et al.|2010| Stello et al.|2011 


Mosser et al.|2012 


)• 



this is our poor theoretical understanding of the underlying 
physical mechanisms for mode driving and damping. 

Using a large set of red giant stars observed by CoRoT, 
Baudin et al. (2011 ) have derived scaling relations in terms 
of mode lifetimes and amplitudes. The se authors have found 
that the scaling relation proposed by Samadi et al. (2007) 
for the mode amplitude significantly departs from the mea- 
sured one. This result was recently confirmed by Huber 



semble asteroseismology by deriving scaling relations that 
relate seismic parameters to a few fundamental stellar pa- 
rameters (e.g. masses, radii, luminosities etc). These ap- 
proaches are now commonly applied to global seismic pa- 
rame ters, such as the cutoff-frequency or peak frequency 



(e.g., 


Miglio et al. 


2009 


Stello et al. 2009 


Kallinger et al. 


2010 


Mosser et al. 


2010 


. However, scaling relation is used 



only infrequently for mode amplitudes. The main reason for 



et al] ( 2011| ), |Stello et al.| ( |2011| ) and |Mosser et al. ( |2Tj12| ) 
with Ke pler ob servations , and is easily understood by not- 
ing that Samadi et al. (2007) established the scaling for 
for main-sequence stars only, and only for mode surface 
velocity. Indeed, those results point out that a dedicated 
theoretical investigation of mode amplitudes in intensity 
for red giants is needed to provide an adequate theoretical 
background. 

Towards the end of their lives, low-mass stars greatly 
expand their envelope to become red giant stars. As a con- 
sequence, the low density of the envelope favours a vig- 
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orous convection such that excitation of solar-like oscilla- 
tions occurs in a medium with very different physical con- 
ditions than encountered in the Sun. This introduces new 
problems about the physical mechanism related to mode 
driving. For instance, the higher the turbulent Mach num- 
ber, the more questionable the assumptions involved in the 
theory (|Go~ldreich fc Keeley] [19771 [Goldreich et al.||1994 



Samadi & Goupil|2001||Chapiin et al.|2005||Belkacem et al 



20101. 

In addition, red giant stars are characterised by high lu- 
minosities and hence have relatively short convective ther- 
mal time-scales at the upper most part of their convective 
envelope. One can therefore expect a stronger departure 
from adiabatic oscillations because the perturbation of en- 
tropy fluctuations related to the oscillations dimensionally 
depends on the ratio L/M (where L is the luminosity and 
M the mass) . Thus, extreme physical conditions in the up- 
permost convective regions of red giants raise new questions 
about the energetic aspects of damped stochastically ex- 
cited oscillations (more precisely mode driving and damp- 
ing). In the present paper, we focus on modelling mode 
driving. We derive scaling relations for red giant stars in 
terms of mode amplitude (in velocity and intensity) and 
compare them with the available CoRoT observations. 

This paper is organised as follows: from a grid of 3D 
hydrodynamical models representative for the upper layers 
of red giant stars, we derive in Sect. |2| theoret ical scaling 
laws for the mo de a mplitudes in velocity (Sect. 2.1) and in 
intensity (Sect. 2.4). These scaling laws are then compared 
in Sect. [3] with seismic data. Finally, Sect. [4] is dedicated to 
conclusions. 



2. Theoretical scaling relations for mode 
amplitudes 

In this section our objective is to compute theoretical scal- 
ing relations of mode amplitudes both in terms of surface 
velocity and intensity. To this end, the mode amplitude will 
be computed with the help of hydrodynamical 3D numeri- 
cal simulations. 

2.1. Surface velocity mode amplitude, v 

The mean-sq uared surface v elocity for each radial mode is 
given by (e.g. |Samadi||2011| and references therein) 



v 2 (is, r) = 



t{v) 9(u) 



2 M(i/,r) ' 



(1) 



where v is the mode frequency, CP the mode excitation rate, 
r the mode life-time (which is equal to the inverse of the 
mode damping rate 77), M the mode mass, and r the radius 
in the atmosphere where the mode velocity is evaluated. 
The mode mass M is defined for radial modes as 



1 



M^rW Jo 



\ti r (v,m)\ 2 dm, (2) 



where £ r is the radial component of the mode eigendis- 
placement. The quantities v, M and £ r are evaluated at 
two relevant layers: 

— the photosphere, i.e. at r — R* where is the stellar 
radius; 



at a layer where spectrographs dedicated to stellar seis- 
mology are the most sensitive. According to |Samadi| 



et al. (2008), for the Sun and solar- type stars, this 
layer is close to the depth where the potassium (K) 
spectral line is formed, that is at the optical depth 
t 500 nm — 0.013. For stars with different spectral type 
this layer may vary, but by an as yet un known manner 
(see the discussion in Samadi et al |2008"| . By default we 



therefore adopt this reference optical depth to be repre- 
sentative for the Doppler velocity measurements for red 
giant stars. This assumption is discussed in Sect. |3.2| 



In Eq. 



Q , y a 

aT|(p008 



CP and M are computed in the manner of 



Samadi et aL| ( 2008 ) using a set of 3D hydrodynamical mod- 
els of the upper layers of sub- and red giant s t ars. H owever, 



this calculation differs from Samadi et al. (2008) in two 
aspects. First, instead of adopting a pure Lorentzian func- 
tion for the eddy-time correlation in the Fouri er domain, 
we introduce, following Belkacem et al. (2010), a cut-off 
frequency derived from the sweeping assumption. Second, 
the 3D models at our disposal have a limited vertical ex- 
tent that results in an under-estimation by up to ~ 10 % of 
the maximum of CP. To take into account the driving that 
occurs at deeper layers we extend the calculation to deeper 
layers using standard ID stellar models (see below). 

The 3D hydro dynamical models were built with the 



CQ 5 BOLD code flFreytag et al. 



(F 
1.2 



2002 



All 3D mode 1 



Wedemeyer et al 
have a solar met al 



2004| Freytag et al.|2012 

abund ance. The chemical mixture is based on |Asplund et al. 
( 2005 ) . The characteristics of these 3D models are given in 
Table [l] All models have a helium abundance of Y = 0.249 
and a metal abundance of Z — 0.0135. The 3D models SI, 
S2, S3, and S7 correspond to red giant stars while S4, S5 
and S6 to sub-giants stars. 

For each 3D model, an associated complete ID model 
(intcrior+surface) is computed in such a way that the outer 



layers are obtained from the 3D model (see Samadi et al. 
for details) while the interior layers are comp uted 



2008 



using the CESAM2K code ( |Morel & Lebretonj [2008| . In 
these ID models, c onvection is treated according to the 
Canuto et al. (1996) local formulation of convection. This 
formulation requires a prescription for the size A of the 
strongest eddies. We assume that A = aH p where H p is 
the pressure scale height and a a parameter adjusted such 
that the interior model matches the associated 3D model 



as detailed in Samadi et al. (20081. The complete models 



(intcrior+surface) are from now on referred to as patched 
models. 

The characteristics of the patched models are given in 
Table [2] We then computed the global acoustic modes 
associated with each of the patc hed models using the 



adiab atic pulsation code ADIPLS ( Christensen-Dalsgaard 
2008). Finally, the mode lifetim es r are evalu a ted u sing 



the meas urements performed by |Baudin et al. 
Sect 



3.11 



(2011 



see 



Our objective is to establish a scaling for the maximum 
Knax hereafter) as a function of stellar pa- 
assuming that t he mo de lifetime r is known. 



of v (Eq. [T 
rameters anc 
As shown by 



Belkacem et al. (2011), the mode lifetime r is 



expected to reach a plateau at a characteristic frequency, 

the maximum of (CP/M) 



As we will see in Sect. 2.2 



also peaks at i/ max . Accordingly, to derive a scaling law 
for V, nax , one needs to determine how t he r atio (CP/M) max 
scales with stellar parameters (see Sect. 2.2). 
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Label log g 



SI 

S2 
S3 
S4 
S5 
S6 
S7 



2.50 
2.50 
2.00 
3.50 
3.50 
3.50 
3.00 



T e ff 

[K] 



4964 
4475 
4551 
4931 
5431 
5885 
5039 



22 
10 
16 
20 
23 
16 
11 



Table 1. Characteristics of the 3D models. T e g is the ef- 
fective temperature, and g the surface gravity. 



Table 2. Characteristics of the ID "patched" models. L is 
the luminosity, M the mass, Av the large separation, and 
v c the acoustic cutoff-frequency. 



Among these parameters, apart from the classical fun- 
damental parameters (luminosity L, mass M, effective tem- 
perature T c g, gravity g, etc), we in addition considered the 
acoustic cut-off fre quency v c and the large frequ ency sepa- 
ration Av (see e.g. |Christensen- Dalsgaard| |1982 ), since the 
former is related to the properties of the surface and the 
latter to the mean density of the star. These parameters 
scale as 



Av = Av G 




(3) 
(4) 



where quantities labelled wit h the symbol refer to solar 
values, v c q — 5100 /xHz (see [Jimenez 2006, and reference s 
therein), and Avq — 134.9 /mz ( jToutain fe Froehlich 1992 1. 
The values of v c and Av associated with each model are 
given in Table[2] 

Finally, we stress that the characteristic frequency v maxi 
at which r reaches a plateau and CP/M is maximum, is re- 
lated to a resonance in the uppermost layers of solar-like 
star s between the thermal time-scale and the modal period 
(see Belkacem et al. 2011 and reference therein). This is 
why it scales as the acoustic cut-off frequency v c in very 
good approximation: 



^max — ^max,0 



y c,0 



(5) 



where t-Wx,© = 3101 /iHz. 



2.2. Scaling relation for (CP/M) 



The maximum of CP is plotted in Fig. [T] (top) as a function 
of the ratio L/M oc T^ s /<?. This dependence with T e g and g 



was already highlighted an d explained by |Stein et al. £2004 ) 



and Samadi et al. (2007) (see also the review by Samadi 



2011 ) , and is nicely confirmed by Fig. [I] (top) . Indeed, CP D 
follows a power law of the form 



J may J m 



L M 



with s = 2.60 ± 0.08 , (6) 



where CP^ ax = (4.2±J;^) x 10 15 J/s. The maximum of CP is 
found to peak at a frequency close to v max . We note also 
that the value of the exponent and the constant CP^ aY in 



Eq^Jp]) are compatible with the results of Samadi et al 
( 2007^ established on the basis of a small set of 3D models 



Label 


M 


Q 


logs 


Teff 


L 


Av 


V c 




[M e ] 






[K] 


[Le] 


[MHz] 


[MHz] 


Ml 


3.74 


0.565 


2.51 


4962 


172.5 


3.43 


63 


M2 


0.98 


0.621 


2.50 


4463 


30.4 


4.77 


67 


M3 


4.20 


0.610 


1.99 


4551 


444 


1.40 


21 


M4 


1.39 


0.636 


3.53 


4927 


5.86 


25.11 


637 


M5 


1.74 


0.596 


3.50 


5392 


11.5 


23.30 


607 


M6 


1.73 


0.576 


3.51 


5856 


15.9 


23.30 


583 


M7 


2.49 


0.615 


3.00 


5040 


39.3 


9.00 


199 



of the surface layers of main-sequence (MS) stars. We thus 
confirm the validity of this relation from MS to red giant 
stars. 

We turn now to the mode mass, M. Because we aim 
to compare theoretical mode velocities with measurements 
made from the ground with spectrographs dedicated to 
stellar seismology, we eva luate M at the optical dept h 
t 500 nm = 0.013 (see Sect. 2.1 and Samadi et al. (2008)). 
For a given model, the mode mass (M) decreases rapidly 
with v, but above a characteristic frequency close to v ma _ x 
it decreases more slowly. Although M does not have a min- 
imum, we found that, as CP, the ratio (CP/M) reaches a max- 
imum close to fc'max, which scales as given by Eq. ^ and 
Eq. ([5]). Therefore, we evaluate M at v = v niax . From now 
on we label this quantity as M max . 

Among the different stellar parameters mentioned in 
Sect. |2.1[ a clear correlation of M max is found with g, 
(L/M), v c or Av. However, the more pronounced correla- 
tion is found with Av. We therefore adopt the scaling with 
Av. The variation of M max with Av is shown in Fig. [T] 
(bottom). M max can be nicely fitted by a power law of the 
form 



M = M° 



Av ' 



with p = 2.1 ± 0.1 , (7) 



where M° lax = (4.5^^3) x 10 21 kg, and Av is given by the 
scaling relation of Eq. Q . 

By using Eq. ^ andEq. Q , the maximum of the ratio 
CP/M then varies according to: 



( W max = (^max/M 



° ) 

max / 



L M Q 



Av 
Av~^ 



(8) 



2.3. Scaling relation for V max 

Equation ([8]) now permits us to proceed by considering the 
scaling lawfor mode amplitudes, in terms of surface veloc- 
ities. The maximum of the mode surface velocity, by using 
Eq. (|8| together with Eq. 0, reads 



Kiax — Vq 



L M Q 



Av 



V 



where r„ 



t M j \Av Qi 

is the characteristic lifetime at v = v„ 



t'o 




(9) 
{ , and 

(10) 



with To a reference mode lifetime whose values are arbitrary 
fixed to the lifetime of the solar radial modes at the peak 
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100 



where 5L(t) is the mode Lagrangian (bolometric) luminos- 
ity perturbation, 5T e g(t) the effective temperature fluctua- 
tion, and 5R*(t) the variation of the stellar radius. 

Since the second term of Eq. ( 11 ) is found negligible in 
front of 5T c ff (i), one obtains the rms bolometric amplitudes 
according to 



SL 
L 



\ 1 eS / , 



(12) 



where the subscript rms denotes the root mean-square. 

We now need a relation between (£T e ff/T eff ) rms (or 
equivalently (SL/L) ims ) and the rms mode velocity V max . 
For convenience we introduce the dimensionless coefficient 
C defined according to 



ST c g 



cir 



= C (SL/L) 





rms 



Vq 



(13) 



where {6L/L)f m 

the solar bolometric mode amplitude (IMichel et al.||2009|), 



2.53 ±0.11 ppm is the maximum of 



TSf = 5777 K the effective temperature of the Sun, and 
^rms = 18-5 i 1-5 cm/s the maximum of the solar mode 
(intrinsic) s urface velocity evalua ted at the photosphere as 
explained in Samadi et al. ( 2010[ ). 

The quantity £ in Eq. ( |13p is defined at an arbitrary 
layer, which is generally the photosphere (i.e. at r = i?*). 
Accordingly, we must evaluate the velocity and hence the 
mode mass M at that layer. This implies the following scal- 
ing for M max : 



TVT — TVT U 

Jn max,t — Ju max,t 



Av 



(14) 



Fig. 1. Top: CPmax as a function of L/M. The triangles are 
associated with the 3D models. The red line is a power law 
of the form (L/M) s with s — 2.6. Bottom: Mode mass 
at v = i/ max (M max ) as a function of the large separation 
Av. The mode masses are evaluated here at the layer cor- 
responding to the optical depth r 50 o n m = 0.013 (see text). 
The triangles are associated with the 3D models. The red 
line is a power law of the form (Av / Avq)~ p with p = 2.1. 



where p. M = 2.0 ± 0.10 , M° lax , = (8.0±|;f) X 10 21 kg and 
Av is given by the scaling relation of Eq. Q. 

Combining Eq. ( 13 ) with Eq. ([9]) gives the scaling for 
the bolometric amplitude 



SL 
~L 



/SL\ Q fvo, 

V L / rms \vSns 



'IM.-.V / L Mq 



T 



Av 

~Av® 



(15) 



frequency, that is To = 3.88 days. Accordingly, we have where vq 
vq = 0.41 m/s. 

It is worthwhile to no te that our scaling relation ( Eq. [9]) 
differs from the result of Kjeldsen & Bedding (2011). This 
is explained by the fact that the postulated relation of 
|Kjeldsen fc Beddin g (2011) for mode amplitudes in velocity 
(their equation 16) does not take the mode masses into ac- 
count, while this is definitively necessary as seen in Eq. (JlJ. 



_ to ( V° m 



2 VM° 



= 0.31 m/s. 



2.4.1. Adiabatic case 



Within the adiabatic approximation, it is possible to relate 
the mode surface velocity to intensity perturbations (e.g., 
Kjeldsen fc Bedding|p95] ) ; this give: 



2.4. Scaling relation for bolometric amplitude 

The instantaneous bolometric mode amplitude is deduced 



CK95 




(16) 



at the phot osphere according to ( Dziemblowski 1977 
Pesnell|[l990l ) 



6L(t) . ST eB (t) , n 8R*(t) 



L 



(11) 



which assumes that the modes are quasi-adiabatic, but 
not only. It supposes that the modes propagate at the 
surface where they are measured. This approximation is 
not valid in the region where the modes are measured 
since in this region they are evanescent. Furthermore, it 
assumes an isothermal atmosphere. A more sophisticated 
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quasi-adiaba tic approach has been proposed by |Severino] 



et al. (2008). The authors went beyond the approximation 
of isothermal atmosphere by taking into account the tem- 
perature gradient as well as the fact that the intensity is 
measured at constant instantaneous optical depth. Both 
effects are taken into account by the method described 
in Sect. |2.4.2[ which in addition considers non-adiabatic 
modes. 

We present in Fig. [2]Ck95 as a function of (L/M). The 
adiabatic coefficient remains almost constant for the type 
of stars investigated here (sub- and red giant stars). This is 
obviously because Ck95 varies as the inverse of the square 
root of T e ff. 



2.4.2. Non-adiabatic case 

We also c omputed ( using the MA D non-adiabatic pulsa- 
tion code (Grigahcene et al. 2005). This code includes the 
time-dependent convect ion (TDC) treatment described in 
Grigahcene et al. (2005). 

This TDC formulation involves a free parameter j3, 
which takes complex values and enters the perturbed energy 
equation. This parameter was introduced to prevent the 
occurrence of non-physical spatial os cillations in the eigen- 
functions (see Grigahcene et al.||2005 for details) . To con- 



L/M 



strain this parameter we used the scaling relation between 
the frequency of the maximum height in the power spec- 
trum (i-'max) and the cut-off frequency (v c ). When scaled 
to the Sun, one can use this scaling to infer ^ max for the 
models we used and the parameter (3 is then adjusted so 
that the plateau (o r depression) of the co mputed damping 
rates coincides (see Belkacem et al.|[2012 ). 

Note also that TDC is a non-local for mulation of con- 
vection an d is based on t he Gabriel ( 1996[) formalism as ex - 
( |2006bp and |Dupret et al] ( |2006a I . 
ocal parameters related to the con- 



plained in Dupret et al. 
In this framework, non- 



vective flux (a) an d the turbulent pressu re (b) are chosen in 
the same way as in Dupret et al. ( 2006b| see their Eqs. (17) 
and (18), see also Dupret et al. ( |2006 c)) so that it fits the 
solar 3D numerical simulation. This calibration results in 
a = 10.4 and b — 2.9 (assuming a mixing-length parameter 
a = 1.62) 

For sub- and red giant stars (L/M > 10 Lq/M & ), the 
non-adiabatic intensity-velocity relation obtained with the 
MAD code can quite well be fitted by a power law of the 
form 



Cnad — Co 



L M Q 
~L^ ~M 



(17) 



where k = 0.25 ± 0.05 and Co = 0.59 ± 0.07. For main- 
sequence stars (L/M < 10i Q /M Q ), Cnad remains almost 
constant (not shown). For the Sun, we find Cnad — 0.95, 
which is close to the value expected by definition for 
the Sun. Therefore, we are then led to multiply Cnad 
by only a factor 1.05 so that, for the Sun, theoretical 
(<5-L/-L) max matches the helioseismic measurements. The re- 
sult is shown in Fig. [2] for the sub- and red giant stars 
(L/M > 10 (Lq/Mq)). The non-adiabatic coefficient in- 
creases rapidly with increasing (L/M) while Ck95 remains 
almost constant. Hence, the higher (L/M), the larger the 
difference between the non-adiabatic and the adiabatic co- 
efficient (Crt95)- 



Fig. 2. Coefficient C (see Eq. 13) as a function of L/M for 
sub- and red giants. The filled circles correspond to the val- 
ues, Cnad) obtained with the MAD non-adiabatic pulsation 
code (see details in the text ). The empty squares corres pond 
to the adiabatic coefficient Kjeldsen & Bedding ( 1995 ) (see 
Eq. 16 ). The red line corresponds to a power law of the form 
Co (L/M) k with fc=0.25. Both intensity-velocity relations C 
have been calibrated so that for the Sun C = 1 (see text). 



3. Comparison with the observations 

We compare in this section theoretical mode amplitudes 
with seismic measurement s ma de from the ground in terms 
of Doppler velocity (Sect. |3.2[) and from space by CoRoT 
in terms of intensity (Sect. p^3| . We recall that computing 
the theoretical mod e amplitudes requires knowledge of T max 
(see Eqs. [9] and 15 1, which is obtained from a set of CoRoT 
targets as explained in Sect. |3.l) 



3.1. The CoRoT data set 



Baudin et al. ( |2011 ) have measured the mode amplitudes 
for 360 CoRoT red giant targets. Among those targets, 
many show very narrow peaks, close to the frequency res- 
olution of the spectrum, while the others have resolved 
peaks. About 65% of those targets have a highest mode 
whose width is sufficiently broad to be fitted with a 
Lorentzian profile. For those targets, the height of the high- 
est mode, -ffmax, and its lifetime T max are thus derived 
from the fit procedure. However, it is not excluded that 
some modes with a width more narrow than the frequency 
resolution may have been fitted with a Lorentzian profile 
because of the low signal-to-noise ratio. To exclude those 
modes, we only considered modes with a width r max = 
l/(7TT max ) broader than twice the frequency resolution of 
the spectra (which is 0.081 /iHz). This subset represents 
about 170 targets for which we have an estimate of the 
mode lifetime (r max ) at the peak frequency. For each tar- 
get of this subset, the maximum of the mode amplitude 
in intensity (^4 max ) was obtained according to the relation 

^4max = \J -ffmax/Tmax- Fin ally, a bolome t ric cor rection was 
applied in the manner of Michel et al. (2009) to convert 
the apparent intensity fluctuation A„ 
amplitude (i5L/L) max . 



into a bolometric 
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3.2. Maximum velocity amplitude (V max J 

The mode amplitude in terms of velocity is given by Eq. Q. 
Calculating V max requires to know the mode life time r max 
at the peak frequency. We used the values of_r max avail- 
able for our set of CoRoT targets (see Sect. 3.1). We also 
determined the ratio L/M as well as Aza The luminosity 



and mass of these targets are unknown. However, |Baudin| 
et al. (2011) have proposed to derive an estimate of the 
ratio L/M using the following scaling: 



L T r 
M X v, 



7/2 



eff 



(18) 



where i/ max is the frequency of the maximum mode height 
H max and T c g is determined fr om photometri c broa d-band 
measurements as explained in Baudin et al. (2011). Note 
that the scaling law of Eq. ( 18 ) assumes that 
as i/ c , which scales as g /y/T e s 

as first established by Stello et al.| (I2T509J), |Hekker et al 



(200! 

relat: 



^max scales 
C oncerning Av, 



2009) and Kallinger et al. (2010), there is a clear scaling 



ation between this quantity and u max . We derived this 
quant i ty her e according to the relation derived by |Mosser| 
et al. (2010) from a large set of CoRoT red giant stars: 



Av = 0.280 



,0.747 



(19) 



Theoretical values of Vmax were compared with the stars 
whose Vmax has been measured so far in Doppler velocity 
from the ground. We considered the different measurements 



published in the literature I 


Frandsen et al. 


2002 


Barban 


et al.||2004| |Bouchy et al.| 


2005| |Carrier et al.| 


2005aE 


Mosser et al.|2005||Arentoft et al. 2008 Kjcldsen et al.|2008 


Mosser et al.||2008l |Teixeira et al.||2009 


Ando et al. 


2010). 



terms of peak amplitudes. In that case they were converted 
into root-mean-square (rms) amplitudes. Furthermore, we 
rescaled all amplitudes into intrinsic (by opposition to ob- 
served) amplitudes. Measured values of Vmax are shown in 
Fig. [3] (top panel) as a function of L/M. We have an esti- 
mate of the ratio L/M for only a few stars while for almost 
all of them we have a seismic measure of u ma]t , which is 
typically more accurate than the determination of the ra- 
tio L/M. Therefore, we also show V max in Fig. [3] (bottom) 
as a function of v max . The theoretical values of K ma3( ob- 
tained for our subset of red giants are found to be close to 
the measurements obtained for the red giant stars observed 
in Doppler velocity from the ground. Note that the con- 
siderable dispersion seen in the theoretical values of V max 
comes from the dispersion in the measured value of r max . 
Furthermore, we point out that the parameters p, s, 7q, 
and Mo, which appear in Eq. ([8]), are mostly determined 
with quite a large error. The errors associated with the pa- 
rameters introduce a bias on the theoretical Vmax, which is 
shown in Fig. [3] by a red vertical bar. As seen in Fig. |3j the 
theoretical Vmax are found, on average, to be about 30% 
lower than the measurements. 

Using sever al 3D simulat i ons o f the surface of main- 
sequence stars, Samadi et al. (2007) have found that Vm ax 
scales as (L/M) sv with sv = 0.7. As seen in Fig. [3| this 
scaling law reproduces the MS stars quite well. When ex- 
trapolated to the red giant domain (L/M > 10 Lq/Mq), 
this scaling law results for Vmax in values very close to our 
present theoretical calculations. 



10.0 
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Fig. 3. Top: Maximum of the mode velocity V max as a func- 
tion of L/M. The filled circles correspond to the MS stars 
observed in Doppler velocity from the ground and the red 
line to the power l aw of the form (L/M) 7 obtained by 
Samadi et al. (2007 1 using 3D models of MS stars. The blue 
squares correspond to t he t heoretical V max derived accord- 
ing to Eq. ([9]) (see Sect.|3.2|). The red square corresponds to 



the median value of the theoretical V max and the associated 
vertical bar corresponds to bias introduced by the l-o~ error 
associated with the parameters p, s, To, and Mo (Eq. [8|. 
Bottom: Same as top as a function of u max . 



so far evaluated a 
= 0.013 (see Sect 



the 



2.1 ■ 

to the optical depth 



reference optical depth r 500 nm : 
We now discuss the sensitivity M 
at which they are computed. To evaluate our sensitivity 
to this choice, we alternatively computed the theoretical 
Viiax a t the photosphere and at an optical depth ten times 
lower than our reference level, that is at r 500 nm = 10 -3 . 
Theoretical Vmax are found to be ~ 30% lower at the 
photosphere and higher by ~ 20% at the optical depth 
7~ 500 nm = 10~ 3 . This result illustrates at which level V max 
is sensitive to the depth where the acoustic modes are sup- 
posed to be measured. This depth is not well known, how- 
ever, but we believe that it should be located between the 
photosphere and our reference optical depth. 
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3.3. Maximum bolometric amplitude ((5L/L) max ) 
3.3.1. Adiabatic case 



We computed (SL/L) max according to Eq. (151 using the 
scaling law given by Eq. (|9| for v and assuming the 
adiabatic coefficient (k95 (Eq. ([16])). Fig. [4] (top) shows 
(SL/L) max as a function of ratio (L/M), where this ra- 
tio is estimated according to Eq. (18). We also plotted the 



mode amplitudes measu red for a small sam ple of CoRoT 
main-sequence stars (see Baudin et al.|201l" and references 
therein). Theoretical (SL/ L) max under- estimates the ampli- 
tudes measured on the CoRoT red giant stars by a factor 
of about 2.5. 



3.3.2. Non-adiabatic case 

We computed (SL/L) max according to Eq. ( |l~5"| ) assum ing 
the non-adiabatic scaling law established in Sect. \2.4~2 (see 
Eq. for C- The result is shown in Fig. [I] (Hottom). 

Using the non-adiabatic coefficient results in an increase 
of the bolometric amplitude by a factor ~ 1.5 compared 
to the calculations based on the adiabatic coefficient. This 
renders the theoretical bolometric amplitude closer to the 
observations. 

We have plotted in Fig. [5] the histogram of the rela- 
tive difference between observed and theoretical (5L I 'L) m „, 
that is, the histogram of the quantity 7 = (A obs — A) /A, 
where A is the theoretical amplitude and A obs the observed 
one. The dispersion seen in the histogram is due both to 
the errors associated with the data and the fact that we 
observe a heterogeneous population of stars with different 
chemical abundance. 

The red horizontal bar shows the bias introduced by the 
1-er errors associated with the determination of the param- 
eters p*, s, To, Mo*, k, and Co as well the measurement 
of (SL/L)f 



and 



(see Eq. (13l). The median of 7 is 



close to 0.8 (the vertical dashed line). This means that the- 
oretical amplitudes remains, on average, ~ 40 % below the 
CoRoT measurements. 



4. Conclusion 

4.1. Theoretical scaling relation for the velocity mode 
amplitude 



We have ex tended the calculations performed by |Samadi| 
et al. (20071 for main-sequence stars to sub- and red giant 
stars. We found that the maximum of the mode excitation 
rate, T max , scales approximately as (L/M) s with s — 2.60± 
0.08. Accordingly, for sub- and red giant stars, theoretical 
Tmax scales in same way as for the main-sequence stars. 

We also found that the mode mass at the peak fre- 
quency, M max , which was evaluated at a reference level in 
the atmosphere, scales as Av~ p where Av cx (M/i? 3 ) 1 / 2 , 
with p = 2.1 ± 0.1. Since (M/R 3 ) represents also the mean 
density, we have that M max scales almost linearly as the 
inverse of the star mean density. This tight relation still 
remains to be understood, however. 

From the scaling laws for M max and Tnia*, we finally 
derived a scaling law for the maximum of the mode velocity, 
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Fig. 4. Top: Maximum of the mode intensity fluctuation 
(SL/L) max as a function of L/M. The filled circles corre- 
spond to th e seismic measurements performed by |Baudin 



et al. 



(2011) on a large number of CoRoT red giant stars 
T?0 targets). We only considered the targets for which 
the mode line widt h is broader than twice the frequency 
resolution (see Sect. 3.1). The empty circles correspond to 
the M S stars observed so far by CoRoT (see Baudin et al. 
2011), and the blue squares are the theoretical (6LJL) max 
computed according to the Kjeldsen fc Bedding| (1995) 



adiabatic coefficient (Eq. (16), see Sect. 3.3.11. The red 



square corresponds to the median value of the theoretical 
(5L/L) max and the associated vertical bar corresponds to 
the bias introduced by the 1-er error associated with the 
parameters p*, s, To, and Mo,*, {&L/L)f ms and vf ms (see 
Eq. ( 15 )). Bottom: Same as top, the theoretical (SL/L) max 
are computed here assuming fo r C th e non-adiabatic scal- 
ing relation established in Sect. 2.4.2 (see also Fig. |2j. The 
red error bar here also accounts for the 1-er error associated 



with the parameters k and Co (see Eq. (17) and Sect. 2.4.2) 



which has the following form: 



Kiax C< (^max) 



1/2 



s/2 



M 
R 3 



p/4 



(20) 



where r max is the mode lifetime at the peak frequency. 



Using CoRoT data, Baudin et al. (2011) have found 



that r max scales approximately as where m = 16.2 ± 
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Fig. 5. Histogram of the relative difference (7) between 
observed and theoretical {$L/ L) (see text). The verti- 
cal dashed line shows the position of the median value. The 
horizontal error bar corresponds to the bias introduced by 
the Tcr error associated with the determination of the pa- 
rameters p*, s, Tq, Mo,*, k, and Co and the measurement 
^ (SL/L)f mB and *,° ms . ' 



2 for the main-sequence and sub-giant stars. Recently, 
Appourchaux et al. ( 2012 ) have found a slope m = 15.5±1.6 
wi th Kepler data , whic h is hence compatible with that 
of Baudin et al. (2011). Such a power law is also sup- 
ported by the theoretical calculations of Belkacem et al. 
(2012) performed for main-sequence, sub- and red giant 
stars. Furthermore, although M max scales better with Av, 
it also scales well as g~ p with p' — 1.66 ± 0.15 (note the 



larger uncertainty for p' compared to p). Accordingly, since 
L/M oc T^ s /g, we can rewrite the scaling for V mSLX (Eq. 20) 
as a function of the star spectroscopic parameters only: 



Vmax oc T 



(2s-m/2) (p'/2-s/2) 



err 



(21) 



Using a set of C0R0T red gi ant stars for which t he mode 
lifetimes have been measured ( |Baudin et al. 2011), we de- 
rived from the scaling law of Fq. theoretical values of 
U m a X . These values were found to be close to the measure- 
ments made from the ground in terms of Doppler velocity 
for red giant stars. However, the Doppler measurements re- 
main on average under-estimated by a about 30%. We dis- 
cuss in Sect. [5]possible reasons for this under-estimation. 



4.2. Theoretical scaling relation for the bolometric mode 
amplitude 



When converted in terms of intensity using the Kjeldsen 
fc Bedding] ( |1995[ ) adiabatic relation, the theoretical ampli- 
tudes un der-estimate the bolo metric mode amplitudes mea- 
sured by Baudin et al. (2011 ) on a set of C0R0T red giant 
stars by a factor about 2.5. Alternatively, we have consid- 
ered the M AD non-adiabatic pulsation code (Grigahcene 
et al. |2005 1 to establish a non-adiabatic relation between 
intensity and velocity. We found that this relation scales as 
(L/M) k with k = 0.25±0.05. We finally established for the 



mode amplitude in intensity the following scaling law: 

s/2+k 



1/2 



M\ p ' /4 



(22) 



where = 2.0 ± 0.1. As for the scaling relation for Vmax, 
the one for (SL/L) max can be rewritten as a function of the 
star spectroscopic parameters only: 

iSL/L) max cx W2+4fc) 9 ^ 2 -^- k ) , (23) 

where p'„ = 1.63 ± 0.15. 

Using the non-adiabatic scaling law for {o~L/L) max re- 
duces the difference between theoretical and measured am- 
plitudes by a factor ~ 1.5. Our analysis hence explains qual- 
itatively the recent results obtained for red giant stars us- 



ing photometric C0R0T and Kepler observations ( Baudin 
et^[20TTl IHuber et al.1|20Tfl |Stello et al.|[20TTj |MosseF 
et al. 112012). Indeed, we stress that theoretical relation ob- 



tained for mode amplitudes in velocity cannot be simply 
extrapolated into photometry because non-adiabatic effects 
dominate the relation between mode amplitude in velocity 
and intensity. 

However, while the non-adiabatic treatmen t imple- 
mented in the MAD code ( Grigahcene et aL]|2005 ) reduces 
the discrepancy with the C0R0T measurements, the latter 
are still underestimated on average by about 40 %. Possible 
reasons for this discrepancy are discussed in Sect. [5j 

5. Discussion 

The mode masses are sensitive to the layer at which they 
are evaluated, which must in principle correspond to the 
height in the atmosphere at which spectrographs dedic ated 



Sect. 3.2 and Samadi et al. [2008) 



tainty associated with the lack of 



to ste llar seis mology are the most s ensitive (See Sect. 2.1 

ir 



However, the uncer- 
knowledge of this layer 



introduces an uncertainty on the co mput ed amplitudes that 
should not exceed ~30% (see Sect. 3.2). 



The discrepancy with the velocity measurements can 
also be attributed to the under-estimation of the mode 
driving. It is not clear which part of the excitation model 
might be incorrect or incomplete. Nevertheless, we believe 
that a possible bias can arise from the way oscillations are 
currently treated in the region where the driving is the 
most efficient (i.e. the uppermost part of the convective 
region). Indeed, in this region the oscillation period, the 
thermal time-scale and the dynamical time-scale are of the 
same order, making the coupling between pulsation and 
con vection stronger and e nergy losses more significant (see 
e.g. Belkacem et al.|[2011 and references therein). We have 
compared non-adiabatic and adiabatic eigenfunctions com- 
puted for the global standard ID model. The non-adiabatic 
eigenfunctions obtained with the MAD pulsation code dif- 
fer from the adiabatic ones only in a small fraction of the 
excitation region. We found a negligible difference between 
excitation rates computed with non-adiabatic eigenfunc- 
tions and those computed with adiabatic eigenfunctions. 
However, we point out that the underlying theory is based 
on a time-dependent version of the mixing-length theory, 
which is well known to be a crude formulation of convection. 
Therefore a more realistic and consistent non-adiabatic ap- 
proach that does not rely on free parameters and that in- 
cludes constraints from 3D hydrodynamical models is re- 
quired. 
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Finally, part of the differences with amplitudes 
(<5L/_L) max measured by CoRoT can be attributed to the 
intensity-velocity relation. Indeed, if we suppose that the 
mode masses are correct, then we must multiply the mode 
excitation rates JVax by a factor <~ 1.5 2 = 2.25 to match 
the velocity measurements. In that case only a differ- 
ence of about 20% with the observed (SL/L) max remains, 
which must then be attributed to the intensity-relation. 
The intensity-relation strongly depends on the way non- 
adiabatic effects are treated, and as mentioned above, the 
current non-adiabatic treatment is based on a crude de- 
scription of the convection and its inter-action with pulsa- 
tion. 
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